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Abstract. We present new methods of automating the construction of in- 
dex pairs, essential ingredients of discrete Conley index theory. These new 
algorithms arc further steps in the direction of automating computer-assisted 
proofs of semi-conjugacies from a map on a manifold to a subshift of finite 
type. We apply these new algorithms to the standard map at different val- 
ues of the perturbative parameter e and obtain rigorous lower bounds for its 
topological entropy for e £ [-7,2], 



1. Introduction 

The literature on computer-assisted proofs in dynamical systems through topo- 
logical methods involves, for the most part, the investigation of the mapping prop- 
erties of certain compact sets which go under different names in different settings: 
they are called index pairs in Conley index theory and windows or h-sets in the 
theory of correctly aligned windows or covering relations [I 104], which are the 
two leading theories applied in topological, computer-assisted proofs in dynamical 
systems. An index pair is not necessarily an h-set and vice- versa, but they both 
satisfy the crucial property that they map across themselves in a way similar to 
the way Smale's horseshoe maps across itself. In other words, they agree in some 
sense with the expanding and contracting directions of the map and thus they can 
be thought of as slight generalizations of Markov partitions. Not surprisingly, after 
several conditions are met in each setting, one is able to prove a semi-conjugacy 
to symbolic dynamics, giving a good description the dynamics of a large set of the 
orbits of a dynamical system, much like one can give a good description of the dy- 
namics of Smale's horseshoe via the easily-describable dynamics of the full 2-shift. 
Given an index pair, an automated procedure was described in [DFT08] to prove a 
semi-conjugacy to a sub-shift of finite type (SFT). 

The first major push towards automated validation and automated computer- 
assisted proofs using topological methods can be found in [DJM05]. By validation 
we mean a mechanism which upholds the claim of a theorem in rigorous way by 
means of finitely many calculations on a computer. Built on those and other results, 
[DFT08] achieved the complete automation of the validation part of a computer- 
assisted proof. In this paper we address the issue of automating the creation of 
index pairs in a computationally efficient way. 
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There are many reasons why automated proof techniques such as [DFT08] come 
at an advantage. For example, for parameter-depending systems, automation allows 
for the rigorous exploration of a system at many parameters. Such applications 
of [DFT08] have already been carried out in [FrolO]. Besides exploration at different 
values of parameters, automated methods may permit us to construct index pairs 
where it would be otherwise impossible to do so by hand. In [DFT08] computer- 
assisted proofs were given of a semi-conjugacy from the classical Henon map to a 
SFT with 247 symbols. These proofs were completely automated, that is, besides 
an input of initial parameters, there was no human intervention from start to finish 
in the computation which proved the semi-conjugacy and bounded the topological 
entropy. 

It should be noted that the complete automation described in the previous para- 
graph is remarkable. To prove a semi-conjugacy using topological methods on a 
computer the following steps must be taken: 

(1) Locate and identify the relevant invariant objects; 

(2) construct a compact set K (refered to above, in our case the index pair) 
whose components map onto each other in a way which is compatible with 
the dynamics of the system; 

(3) prove that some of the dynamics of the system can be described via the 
mapping properties between the components of K. 

Each of these steps presents its own challenges and difficulties. In [DFT08, §3.2] 
step (3) was completely solved, in the sense that given any compact set (index pair) 
K which satisfies the properties required by (2), a completely automated routine 
was given to prove a semi-conjugacy. Steps (1) and (2), the creation of the index 
pair, were also fully automated in [DFT08, §3.1] (and [FrolO]), largely due to the 
properties of "well-behaved" maps like the Henon map: it has a low-dimensional 
attractor and all its periodic orbits exhibit hyperbolic-like behavior from a compu- 
tational perspective. In other words, in such cases there is a dominating invariant 
set, namely the attractor whose periodic orbits all seem hyperbolic and are easily 
localizable. Thus, it is of interest to consider the case where we do not have such 
underlying structure. Quoting directly from [DFTOX, §5]: "... further analysis and 
optimization of the procedure described in section 3. 1 for locating a region of interest 
should lead to even stronger results." 

This article is an effort to extend the methods and techniques presented in [DFT08] 
to more general situations with the aim of dealing with systems which do not ex- 
hibit Henon-like behavior. As an illustration of the approach we propose, we apply 
our method to the standard map 

f e : (x, y) h-> ( x + y + — sin(27rx) mod 1 , y + — sin(27rx) mod 1 ) , 

for (x, y) £ S 1 x S* 1 = T 1 , at different values of the perturbation parameter e. By 
treating the perturbation parameter as an interval, we are also able to give a de- 
scription of the orbits of the map in terms of symbolic dynamics for all perturbation 
values inside the given interval. We chose the standard map for a couple of reasons. 
First, it is a very well-known system with rich dynamics. Second, in contrast with 
previous results for Henon-like systems (e.g. [DFT08, FrolO]), the standard map 
is a volume-preserving map which exhibits both hyperbolic and elliptic behavior, 
which make the task of automating the construction of index pairs considerably 
more difficult. Finally, we generalize the methods of [DFT08] to work on manifolds 
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with nontrivial topologies and not just fi", and our application to the standard 
map exemplifies this. 

In the present constructions, we benefit much from some a-priori information 
of the map. The distinction between a hyperbolic or elliptic periodic orbit, for 
example, is a crucial one. Symmetries of the map, the knowledge of the precise 
location of any heteroclinic or homoclinic orbits are also valuable data from which 
to begin. In our view, a little previous knowledge goes a long way, as we exploit 
this a-priori information to make the constructions highly efficient while keeping 
the entire process highly automated. 

We mention that besides [DFT08] there have been other recent efforts towards 
the automation of computer-assisted proofs in dynamical systems. In particu- 
lar, [MJ09] is a collection of set-oriented algorithms for the creation of index 
pairs adapted for volume-preserving maps. In addition, [BCG07, BCGH08] use 
the method of covering relations after reducing the problem of finding suitable sets 
and maps between sets to a problem of global optimization theory. 

In contrast to other Conley index methods [DFT08, MJ09] which are top-down, 
in that they start with a large portion of the phase space and try to narrow down 
to the invariant objects, our method is bottom-up, in that we are guided by the 
dynamics of the system as we "grow out" the index pairs given good numerical 
approximations of the invariant sets (see Algorithms 2 and 3). The goal of our 
bottom-up approach is to compute a minimal number of images of a discretized 
version of the map (see section 2.3), which is often the most expensive part of 
computations. Our approach thus achieves high efficiency in the automation of the 
construction of an index pair. 

The present work has been motivated by the challenges brought by the standard 
map to the task of automating the construction of index pairs and by what we 
expect will be obstacles for construction of index pairs in future applications. The 
techniques of this paper along with those of [DFT08] and the references therein 
provide a solid set of tools for automating computer assisted proofs in the paradigm 
of planar maps with positive entropy through discrete Conley index theory. These 
methods are technically not restricted to dimension 2, but remain mostly untested 
in higher dimensional examples with higher dimensional unstable bundles. To the 
authors' knowledge, there is an example in [DJM04] with a two dimensional unstable 
bundle, but no others in the existing computational Conley index literature. Having 
focused on making the construction of index pairs highly efficient in this paper, it 
remains to apply them to such higher dimensional examples, which we plan to do 
in future work. 

To the authors' knowledge, the rigorous bounds presented here on the topological 
entropy for the standard map are the only ones of their kind in the literature. We 
could only find [YT09] for computational lower bounds of the topological entropy 
for the standard map at various parameter values, but the results there are non- 
rigorous. Thus it would be interesting to see how other topological methods (for 
example [/ 104, NBGM08]) perform for the same values we have examined in this 
paper. 

This paper is organized as follows: in section 2 we review the necessary back- 
ground for the subsequent sections. Section 3 details the data structures we use in 
our algorithms as well as presents the two new algorithms which we consider as the 
main contribution of this work, Algorithms 2 and 3. We apply these algorithms in 
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section 4 to the standard map to give lower bounds for the topological entropy at 
different values of the perturbative parameter. Theorem 4.1 gives lower bounds of 
the topological entropy of the standard map as given by topological entropy of sub- 
shifts of finite type, for all f £ with e £ [0.7, 2.0]. Theorem 4.2 gives a positive lower 
bound for the topological entropy for e = \. Theorems 4.3 and 4.4 give a positive 
lower bound for the topological entropy for e = 2 by finding connecting orbits to 
hyperbolic periodic orbits of higher period. We conclude with some comments and 
remarks on the implementation in Section 4.4. 

Acknowledgements. We would like to thank Jay James for many helpful discus- 
sions during the course of this work, as well as for his help with the implementation 
of the methods of [CFdlL05] for the standard map, as well as Sarah Day for helpful 
comments on an early draft of the paper. We also thank the Center for Applied 
Mathematics at Cornell University for allowing us access to their computational 
facilities, where the computations in this paper were carried out. We would like to 
thank an anonymous referee for helping us make the paper clearer. 

2. Background 

We review in this section the necessary facts about topological entropy, symbolic 
dynamics, and the discrete Conley index which will be relevant in the later sections. 
Although the exposition will not be in depth, we encourage the interested reader to 
consult section 2 of [DFT08] for a more detailed exposition. For deeper treatment 
of the discrete Conley index, see [MM02]. 

2.1. Topological Entropy and Symbolic Dynamics. Let / : X — > X be a 

map. The topological entropy of / is the quantity h(f) £ K U {oo} which is a good 
measure of the complexity of /. A map with positive topological entropy usually 
exhibits chaotic behavior in many of its trajectories. 

Definition 2.1. Let / : X — > X be a continuous map. A set W C X is called 
(n, e)-separated for / if for any two different points x,y € W, \\f k (x) — f k (y)\\ > £ 
for some k with < k < n. Let s(n,e) be the maximum cardinality of any (n,e)- 
separated set. Then 

h(f\ v \ogs{n,e) 
n{jj = suphmsup 

is the topological entropy of /. 

Although it is always well-defined, it is usually impossible to compute the topo- 
logical entropy directly from the definition. The set of systems for which it is 
computable is rather small, but it contains all systems from symbolic dynamics. 

Let = {0, . . . ,N — 1} Z be the set of all bi-infinite sequences on N symbols. 
It is well-known that E^r is a complete metric space. Define the full N- shift a : 
Ejv — > E^r to be the map acting on Sjv by (<r(x))i — ajj+i. Given an N x N matrix 
A with Aij £ {0, 1} we can define a sub-shift of finite type defined on E^i C E_/v, 
where x £ E^ if and only if for x — (. . . , Xj, a^i+i, ■ ■ ■ ), A x . Xi+1 — 1 for all i. In 
other words, 'Ea consists of all sequences in Ejv with transitions of a allowed by A. 
Viewed as a graph with N vertices and with an edge from vertex i to vertex j if 
and only if Aij — 1, can be identified with the set of all infinite paths in such 
graph. 

The topological entropy of a sub-shift of finite type is computed through the 
largest eigenvalue of its transition matrix A. 
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Theorem 2.2. Let a : Y*a — > £U be a subshift of finite type. Then its topological 
entropy is 

h(a) = logsp(A), 
where sp(A) denotes the spectral radius of A. 

In order to study a map / of high complexity, it is sometimes possible to study a 
subsystem of it through symbolic dynamics. This is done through a semi-conjugacy. 

Definition 2.3. Let / : X — > X and g : Y — > Y be continuous maps. A semi- 
conjugacy from / to g is a continuous surjection h : X — > Y with ho f = g oh. We 
say / is semi-conjugate to g if there exists a semi-conjugacy. 

Note that through a semi-conjugacy h : X —> Y, information about g acting on 
Y gives information about / acting on x. In particular, the topological entropy of 
g bounds from below the topological entropy of /. 

Theorem 2.4. Let f be semi-conjugate to g. Then h(f) > h{g). 

Corollary 2.5. Let f be semi-conjugate to a : Y,a — > ^a, for some N x N matrix 
A. Then 

h(f) > sp(A). 

2.2. The Discrete Conley Index. Let / : M — > M be a continuous map, where 
M is a smooth, orientable manifold. 

Definition 2.6. A compact set K C M is an isolating neighborhood if 

lnv(K,f) clnt(K), 

where Inv(i(f, /) denotes the maximal invariant set of K and lnt(K) denotes the 
interior of K. S is an isolated invariant set if S = Inv(X, /) for some isolating 
neighborhood K. 

Definition 2.7 ( [RS88]). Let S be an isolated invariant set for /. Then P = 
(Pi, Pq) is an index pair for S if 

(1) Pi\Po is an isolating neighborhood for S. 

(2) The induced map 

f M _j f(x) if*, f(x)eP 1 /P , 

JpKX) - I [Pq] otherwise, 
defined on the pointed space (Pi/Po, [Po]) is continuous. 

Definition 2.8 ( [WI170]). Let G, H be abelian groups and ip : G — > G, ^ : P — > P 
homomorphisms. Then 99 and t/> are s/iz/i equivalent if there exist homomorphisms 
r : G — ?> P and s : P — > G and a constant fc 6 N such that 

r o ip — i{j o r, s o ijj — ip o s, r o s — ip k , and s o r — ip k . 

Shift equivalence defines an equivalence relation, and we denote by the class 
of all homomorphisms which are shift equivalent to ip. 

Definition 2.9 ( [FROO]). Let P = (Pi,Po) be an index pair for an isolated 
invariant set S = Inv(Pi\P ,/) and let fp* : P*(Pi,P ;Z) P*(Pi,P ;Z) be the 
map induced by /p on the relative homology groups P*(Pi, Po; Z). The Conley 
index of 5 is the shift equivalence class [/p*] s of f Paf . 
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One can think of two maps /p* and gp^ as being in the same shift equivalence 
class if and only if they have the same assymptotic behavior. Since an index pair for 
an isolated invariant set is not unique, the Conley index of and isolated invariant 
set does not depend on the choice of index pair. There are two elementary yet 
important results which link the Conley index of an isolated invariant set with the 
dynamics of /. The first one is the so-called Wazewski property: 

Theorem 2.10. If [/p,] s ^ [0] s , then S ^0. 

Since similar assymptotic behavior relates two different maps in the same shift 
equivalence class, it is sufficient then to have a map /p* not be nilpotent in order 
to have a non-empty isolated invariant set. In practice, non-nilpotency can be 
verified by taking iterates of a representative of [/p*] s until non-nilpotent behavior 
is detected. We also have an elementary detection method of periodic points. 

Theorem 2.11 (Lefshetz Fixed Point Theorem). Let /* be a representative of 
[fp*]s with maps fk ■ i?fc(Pi, Po;1t) — > Hk(Pi, Pq',%) which are represented by ma- 
trices. Then if 

Hf*) = Et- 1 )* ^/fe^o, 

fe>0 

then f has a fixed point. Moreover, if A(/™) ^ 0, then f has a periodic point of 
period n. 

Since traces are preserved under shift equivalence, A(/*) is independent of the 
representative of [/p*] s , and we may even denote it as A([/p»] s ). 

Corollary 2.12. Let K C M be the finite union of disjoint, compact sets K\, . . . , K„ 
and let S — Inv(K, f) . Let S' = Inv(K\ , fK m ° ' ■ ■ ° fx x ) C S where fxi denotes 
the restriction of f to K{. If 

A([JW m °---°M]a)^0 

then fK m o - - • o /jfj contains a fixed point in S' which corresponds to a periodic 
orbit of f which travels through K\ , . . . , K m in such order. 

This is a strong and useful tool for proving the existence of periodic orbits. How- 
ever, we may have A([//f m o • ■ • o /r-Js) = while [fic m °' " " fKi)* is n °t nilpotent 
and thus have an invariant set which may be of interest since it behaves like a 
periodic orbit. So we have an analogous result based on the Wazewski property. 

Corollary 2.13. Let K C M be the finite union of disjoint, compact sets K%, . . . , K„ 
and let S = Inv{K, f) . Let S' = Inv(K\ , fx m o • • ■ o fx x ) c S where fxi denotes 
the restriction of f to Ki. If 

[(/jf m °-"°/*x).]a?M0]., 

then S' is nonempty. Moreover, there is a point in S whose trajectory visits the 
sets Ki in such order. 

Both corollaries can be useful in different settings when used to to prove sym- 
bolic dynamics through computer-assistance. When existence of periodic points 
is of primary concern, one can use Corollary 2.12. When existence of trajectories 
which shadow a certain prescribed path given by a symbol sequence, but that do 
not necessarily correspond to periodic orbits, Corollary 2.13 is the better tool. This 
may occur when there are higher-dimensional invariant sets to which the restriction 
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of the dynamics is quasi-periodic. In either case, what makes the implementation 
possible is the ease of computability of the traces of the induced maps on homol- 
ogy for Corollary 2.12 and a sufficient condition for non-nilpotency in the case of 
Corollary 2.13. For more details on the implementation, see [DFT08]. 

2.3. Combinatorial Structures. All concepts from the discrete Conley index 
theory from the previous section have analogous definitions in a combinatorial set- 
ting for which computer algorithms can be written. 

Definition 2.14. A multivalued map F : X =t X is a map from X to its power set, 
that is, F(x) C X. If for some continuous, single- valued map / we have f(x) € F(x) 
and F is acyclic, then F is an enclosure of /. 

The reason multivalued maps and enclosures are used in our computations is 
that if they are done properly, they give rigorous results. Moreover, if F is an 
enclosure of / and (Pq, Pi) is an index pair for F (as we will define below), then it 
is an index pair for /. It also follows that if we can compute the Conley index of 
F and process the information encoded in it, we may obtain information about the 
dynamics of /. 

We begin by setting up a grid Q on M, which is a compact subset of the n 
dimentional manifold M composed of finitely many elements Bi- Each element 
is a cubical complex, hence a compact set, and it is essentially an element of a 
finite partition of a compact subset of M. In practice, all elements of the grid 
are rectangles represented as products of intervals (viewed in some nice coordinate 
chart), that is, for Bi € Q, Bi = Ylk=i [ x fe) Vk\- We refer to each element of Q as a 
box and each box is defined by its center and radius, i.e., Bi — (c,, r*j), where Ci and 
Ti are n-vectors with entries corresponding to the center and radius, respectively, 
in each coordinate direction. In practice we have Ti the same for all boxes in Q but 
this is in no way necessary and there may be systems for which variable radii for 
boxes provide a significant advantage. For a collection of boxes /C C G, we denote 
by \K.\ its topological realization, that is, its corresponding subset of M. From now 
on we will use caligraphy capital letters to denote collections of boxes in Q and by 
regular capital letters we will denote their topological realization, e.g., \B%\ = Bi. 

The way we create a grid is as follows. We begin with with one big box B such 
that \B\ encloses the area we wish to study. Then we subdivide B d times in each 
coordinate direction in order to increase the resolution at which the dynamics are 
studied. The integer d will be refered to as the depth. Thus working at depth d gives 
us a maximum of 2 dn (n = dim M) boxes with which to work of each coordinate of 
size 2~ d relative to the original size of the box B. 

Definition 2.15. A combinatorial enclosure of / is a multivalued map J- : Q =t G 
defined by 

T(B) = {B 1 eG: \B'\ HF(B)^0}, 
where F is an enclosure of /. 

In practice, combinatorial enclosures are created as follows. One begins with 
B G Q and defines F(x), x G B = \B\, as the image of B using a rigorous enclosure 
for the map /. Rigorous enclosures are obtained by keeping track of the error terms 
in the computations of the image of a box and making sure the true image f(B) is 
contained in |.F(<B)|. In this paper, all rigorous enclosures will be obtained with the 
use of interval arithmetic using Intlab [Ru i9]. Note that | J*| becomes an enclosure 
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of /. T : Q =4 Q can also be represented as a matrix T with entries in {0, 1} with 
71 j — 1 if and only if Bj G J-(Bi). Moreover T can be viewed as a directed graph 
with vertices corresponding to individual boxes in Q and edges going from box i to 
box j if and only if 7\j = 1. 

Definition 2.16. A combinatorial trajectory of a combinatorial enclosure J- through 
B G G is a bi-infinite sequence 7g = (. . . ,B-i,Bq,Bi, ...) with Bo = Bn G 
and S„ + i G J-(B n ) for all n G Z. 

The definitions which follow are by now standard in the computational Conley 
index literature. We will state definitions and refer the reader to [DFT08] for the 
algorithms which construct the objects defined. 

Definition 2.17. The combinatorial invariant set in J\f C Q for a combinatorial 
enclosure J- is 

Inv(jV, J) = {6 € 5 : there exits a trajectory C A/"}. 

Definition 2.18. The combinatorial neighborhood or one-box beighborhood of S C 
is 

o(6) = {S'Gg : |6'|n|6|^0}. 

Definition 2.19. If 

o(Inv(A/", J 7 )) C N 
then A/" C 5 is a combinatorial isolating neighborhood for J 7 . 

A procedure for creating a combinatorial isolating neighborhood is discussed in 
Section 3 and an given by Algorithm 1. Once we have a combinatorial isolating 
neighborhood, it is possible to define and create a combinatorial index pair. 

Definition 2.20. A pair V = (Vi,Vo) C Q is a combinatorial index pair for the 
combinatorial enclosure J- if its topological realization Pj = \Vi\ is an index pair 
for any map / for which J- is an enclosure. 

One of the main goals of this paper is a procedure to efficiently compute com- 
binatorial index pairs. This is the content of the next section. We now have made 
all definitions necessary to define the Conley index. Computing the induced map 
on homology at the combinatorial level is no trivial task. We refer the enthusiastic 
reader to [KMM04] for a very thourough exposition on computing induced maps 
on homology. In practice, we use the computational package homcubes, part of 
the computational package CHomP [cho], which computes the necessary maps on 
homology to define the Conley index. 

3. Automated Symbolic Dynamics 

The results of [DFT08] show that it is possible to create an automated procedure 
to rigorously prove a semi-conjugacy from a map f : M —¥ M to a subshift E^. 
The procedure detailed in [DFT08] can be summarized as follows: 

(i) Start with a rectangle in M" and obtain a grid of boxes of constant radius by 
partitioning it in each coordinate direction d times for some fixed d. 

(ii) Compute the combinatorial enclosure J- of / in a neighborhood of the area 
of interest, given by the interval arithmetic images of boxes in the grid. For 
some fixed k, create a collection C Q of boxes which are periodic of period 
n < k under T by finding non-zero entries of the diagonal of T n . 
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(iii) Using shortest path algorithms, and denoting by T>ij any shortest path be- 
tween Bi G Vk and Bj G Vu in Q (if there is no such path, T>^ — 0), create a 

collection A C Q by A = Vu U (U»yy ^ij^j ■ Using the appropriate algorithms, 
create a combinatorial isolating neighborhood for A, a combinatorial index 
pair, and compute its Conley index [/p*] s . 

(iv) Working with a representative of [fp*] s as matrices (one for each level of 
homology) over the integers, we filter out all the nilpotent behavior and find 
a smaller representative which exhibits only recurrent behavior. 

(v) We study [/p*] s by its action on the generators of H\{P\, Pq;Z) grouped by 
disjoint components of Pi\Pq, and using Corollary 2.13 we perform a finite 
number of calculations to prove a semiconjugacy from / to E^, where A is a 
to x to matrix, and to is a number less than or equal to the number of disjoint 
components of P\\Pq. 

This method was applied to the Henon map and a semi-conjugacy to a subshift 
on 247 symbols was proved, which gave a rigorous lower bound on the topological 
entropy. 

We make a few remarks about the approach. First, although all periodic orbits 
of the Henon map exhibit hyperbolic behavior, this is not true in all systems. In 
Hamiltonian systems, for example, one expects roughly half of the periodic orbits 
to be isolated invariant sets. Thus looking at the diagonal of T™ is not enough 
to capture orbits which will give us isolated invariant sets and a nontrivial Conley 
index. Another issue which arises in practice is the computational complexity of 
the algorithms to prove a semiconjugacy to a subshift. The complexity increases 
exponentially with the dimension of the system and if we have any hope of getting 
results in systems of dimension higher than 2, we must perform the computations 
as efficiently as possible. 

With this in mind, we view the entire process of proving a semi-conjugacy as the 
composition of two major steps: 

(1) The gathering of recurrent, isolated invariant sets. This can be done 
through a combination of non-rigorous numerical methods, graph algo- 
rithms, set-oriented methods, linear algebra operations, et cetera. 

(2) A proof that this behavior in fact exists through a semi-conjugacy to a 
subshift. 

We point out that the second step is solved in [DFT08]. That is, given as input 
an index pair (Pq,Pi) and the combinatorial enclosure T of / used to create it, 
proving semi-conjugacy from f to a subshift is completely automated. The first step 
is also largely dealt with, but the approach is not as general as it can be, as it was 
devised to deal with maps like the Henon map where isolation is easy. We wish to 
consider cases where one cannot simply compute the map on all boxes in the area 
of interest, either because the dimension of the attractor is too high or because one 
needs a very high box resolution to achieve isolation. Whereas [DFT08] focused 
on (2), here we focus on (1) and on how to produce index pairs more efficiently. 
In settings where isolation is difficult, our methods in this paper to efficiently deal 
with (1) and the solution in [DFT08] of (2) constitute a better method to prove 
semi-conjugacies to a SFT. 

We provide an efficient method of computing an index pair in step (1) when 
one knows roughly what to look for. In essence, we believe that a small amount 
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of a-priori knowledge goes a long way. In other words, we can simplify and make 
computations much more efficient if we know something about dynamics before- 
hand, like previous knowledge of the location of periodic, heteroclinic, homoclinic 
orbits, special symmetries of the system, et cetera. We will illustrate this approach 
with examples in Section 4. 

More precisely, we take as input a list of periodic points and pairs of points which 
might have connecting orbits between them, and produce an index pair containing 
all of the points and connections, if possible. For example, if one had numerical 
approximations of two fixed points and conjectured that there was a connecting 
orbit between them, our algorithms could potentially prove the existence of the 
connection. 

We break the task of finding an index pair into two steps. The first is to approx- 
imate numerically the invariant objects whose existence we want to prove, and the 
second is to cover such an approximation with boxes and add boxes until the index 
pair conditions are satisfied. Before describing these steps, however, we introduce 
our underlying data structures. 



3.1. Data Structures. The algorithms we present rely on some basic data struc- 
tures to encode the topology of the phase space and the behavior of the map. We 
require two routines to keep track of the topology: finding the box corresponding to 
a particular point in the phase space (to compute the multivalued map on boxes), 
and determining which boxes are adjacent (to determine isolation). We then of 
course need a method of storing the multivalued map itself. 

To accomplish all of this, we consider a subset S of our grid and enumerate the 
boxes in S, giving each a unique integer. We then store the position information 
in a binary search tree, which we simply call the tree, so that the index i for 
the box hi covering a given point x in the phase space (i.e., x G \bi\) may be 
computed efficiently. We use the GAIO implementation for our basic tree data 
structure [DFJ01] only, and make use of some enhancements, discussed below. To 
encode the topology, we store the adjacency information in a (symmetric) binary 
matrix Adj, where Adjij = Adj ji — 1 if and only if boxes hi and bj are adjacent. 
More precisely, Adj ij = 1 if and only if \bi\ H | bj | ^ . We will call this matrix 
the adjacency matrix. This is a crucial component when working on spaces with 
nontrivial topology such as the torus, as we will see in section 4. Finally, we store 
the multivalued map as a transition matrix P such that Pjt — 1 if and only if 

bj e F{h). 

Below are the operations of the tree data structure: 

(i) S — boxnums() - Return a list of all box numbers in the tree 

(ii) S — find(C) - Given a set C C M of points, return the indices S of boxes 
covering any i£C 

(iii) insert (C) - Insert boxes into the tree covering points C C M (if they do 
not already exist) 

(iv) delete(S') - Remove boxes with indices in S from the tree 

(v) subdivide () - For each box b in the tree, divide each of its coordinate 
directions in half creating 2™ smaller boxes, thus increasing the depth of the 
tree by 1 
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Although the GAIO tree data structure is in many ways well-suited to our task, 
it has a few simple drawbacks which significantly restrict our computations. Specif- 
ically, the operations insert (•) and delete(-) scramble the box numbers of the tree as 
a side-effect. Thus, if we want to keep track of our current box set while inserting 
or deleting boxes, we have to spend extra time computing the new box numbers. 
Without knowledge of how the renumbering is done, this would take O(nlogn) 
time, where n is the number of boxes, since for each box we have to search the 
tree to find its new number. The best one could hope for would be O(logn) with 
careful bookkepping (one still needs logarithmic time to locate the place in the 
tree for insertion or deletion). We managed to compute the new box numbers in 
0(n) time by observing that the numbering is given by a deterministic depth- first 
search (DFS) traversal of the tree from the root box, and backsolving the modified 
numbers accordingly. 

We make use of several other subroutines as well: 

• P — transition_matrix(tree,/[,S']) - Compute the multivalued map J- 
using / (optionally only for boxes in S) and return it as a matrix 

• Adj — adjacency matrix(tree) - Using the tree.find(. . . ) function, com- 
pute the adjacency matrix which describes the topology of the boxes in the 
tree 

• oS — tree_onebox(S , ,^4c?j) - The indices of boxes in a one-box neighbor- 
hood of S in the tree 

• tree = insert onebox(tree, 5*) - Using tree.insert(. . . ), add boxes to the 
tree that would neighbor a box in S were they already in the tree 

• tree = insert_image(tree,/,5) - Compute the images of each box in S 
under / using interval arithmetic 

• S — grow _isolating(S',P, Adj) - Compute an isolating neighborhood of S 
in the tree. This procedure is called grow dsolating_neighborhood in [DFT08]; 
for completeness we restate it here as Algorithm 1. 

• S — invariant_set(A,P) - Compute the maximal invariant set containing 
N according to the multivalued map P 



Algorithm 1 grow -isolating: Growing an isolating neighborhood 

Input: S,P,Adj 
loop 

/ = invariant-set^, P) 

ol = tree_onebox(/,Ad?) 

if ol C S return I 

S = oI 
end loop 
Output: S 



Note the difference between insert_onebox(tree,S') and tree.onebox^Acij); the 
former adds all boxes to the tree that touch S, but the latter find boxes that touch 
which are already in the tree. 

3.2. Numerical Approximations. We assume here that one has numerical ap- 
proximations of points which are part of hyperbolic, invariant sets such as periodic 
points and homoclinic points. The goal then is to find the desired connecting orbits 
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between specified points. More formally, given a finite set X = {(xi, yi)}i cMx M 
of pairs of points in the phase space, we want to find a small set of boxes S at depth 
d such that [JX C | *S" | and for every pair pi, if Xi € \b x \ and yi €E \b y \, then there is 
a path from b x to b y according to the multivalued map J- . In other words, we wish 
to find connections between each Xi and yi at depth d. Intuitively, we think of the 
set X as the "skeleton" of the invariant behavior of interest. For example, if we 
wanted to connect a period 2 orbit a, b to a fixed point c and back, we could have 
X = {(a, c), (6, c), (c, b), (c, a)}. We assume of course that the desired connections 
exist. 

Typically the most expensive part of the calculations, especially in applications 
involving rigorous enclosures, is the box image calculations. Thus we wish to find 
an algorithm which computes as few images as possible, but also has a reasonable 
running time in terms of the other parameters. 

The first algorithm one might try is to go to depth d, add all boxes in the 
general area of the invariant objects, and then use shortest path algorithms to find 
the connections. Unfortunately this requires computing 0(2 dn ) image calculations, 
where n = dim(M). We can improve this bound by instead doing a breadth-first 
search at depth d, that is, starting at each Xi we compute images until we reach 
yi in a breadth-first manner. This gives roughly 0(|A|6^), where b is the average 
image size of a box, which will typically depend exponentially on n, and £ is the 
average connection length. This is no longer exponential in d, and \X\ is typically 
relatively small, but there are still many image calculations as b is often very large. 

Fortunately we can achieve only 0(d\X\£) image calculations using Algorithm 2, 
which is essentially an recursive version of the first algorithm: we compute the 
connections at a low depth, then subdivide to get to the next depth and grow one- 
box neighborhoods until the connections are found again, and repeat until we reach 
the final depth. 



Algorithm 2 The connection insertion algorithm 

Input: /, X = {(xi,Vi)}i C M x M 
for depth = d sta rt to d cnd do 

tree.insert(Ui{a;i, yi}) 

loop 

for all i do 

Pi = shortest path from Xi to yi in P (or if no path) 
end for 

if Vi pi 7^ 0, break loop 

tree = insert_oiiebox(tree,tree.boxnums()) 
P = transition_matrix(tree,/) 
end loop 

tree.delete(tree.boxnums() \ UiPi) 
tree.subdivideQ 
P = transition_matrix(tree,/) 
end for 



By a strict reading of Algorithm 2, we might end up computing many box images 
repeatedly from the transition_matrix(-) call in the inner loop. To avoid this, we 
simply cache box image calculations: each time a box image calculation is called 
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for, we first check to see whether we have computed it already; if so we look it 
up, and if not we compute it and store it. With caching, we achieve the 0(<i|A|i?) 
bound on image calculations. 

Note that a box connection found at a depth d may only correspond to a eo2~ d - 
chain rather than an actual orbit, where eo is the length of a box diagonal at depth 
0. Typically we rely on shortest path algorithms to give us the most plausible orbits, 
but there are cases, such as spiralling behavior and other 'roundabout' connections, 
where the shortest paths will not be true orbits. For example, even the identity 
map has a path from any box to another according to the transition matrix. In 
all of these cases it may be necessary to use other methods for approximating 
or connecting orbits, such as computing the stable and unstable manifolds of the 
hyperbolic invariant sets to high accuracy and computing heteroclinic intersections. 

Roughly speaking, a long connection between p and q will have many of these 
eo2~ d -chains from p to q to compete with, the vast majority of which do not cor- 
respond to actual trajectories. In fact, in many cases the chains will actually be 
shorter than the true connection. For example, in a simple slow rotation about the 
origin (r, 6) > (r,0 + 2n/k), every point is either period A; or 1 (the origin), but 
one would find closed chains of length much lower than k when sufficiently close 
to the origin. The result of this observation is that longer connections will require 
more iterations of the insert_onebox(tree,tree.boxnums()) call, since the probabil- 
ity of picking a eo2~ d -chain which corresponds to an actual orbit becomes quite 
small as the connection length grows. Consequently, longer connections will take 
significantly longer to compute than shorter ones using Algorithm 2. 

One way around this issue is to make use of the duality between finding con- 
nections between periodic orbits of low period and finding periodic orbits of higher 
period. On the one hand, by finding enough connecting orbits (enough to cap- 
ture horseshoe dynamics) between periodic orbits of low period, one finds infinitely 
many periodic orbits of higher periods. On the other hand, it is often the case 
that by finding enough periodic orbits of all periods (up to some high period) at a 
given depth, most of the connecting orbits which live close to them should also be 
captured. Since dealing with long connections can be problematic, it is simpler to 
deal with specific, high period, periodic orbits if they are easily obtainable. 

A second way to compute long connections is to exploit some a priori knowledge 
of the map. We apply both of these techniques in section 4. 

3.3. Growing an Index Pair. Given a small starting set S of boxes corresponding 
to numerical guesses of hyperbolic invariant sets, we now wish to create an index 
pair from S. To do this, we compute an isolating neighborhood N of S. A first 
approach might be to use the grow -isolating algorithm, Algorithm 1, on the whole 
of M, but as we are concerned with the setting where isolation is difficult, this will 
simply fail to find the desired structures. Even if we manage to cover only what we 
are interested in, this top-down approach requires computing many images; as with 
Algorithm 2, we wish to minimize the number of image calculations while keeping 
a reasonable running time. Specifically, we would like an algorithm that computes 
the images of N and no other images, which would be optimal in our setting. We 
will see that Algorithm 3 does precisely that. 

Note that the loop invariant (and thus correctness) of this algorithm is highly 
sensitive to the order of the steps. As before, we cache box image computations. 
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Algorithm 3 Growing an isolating neighborhood by inserting boxes 

Input: tree,/ 

B = tree.boxnumsQ 

S = B 

repeat 

P = transition_matrix(tree,/,5') 
Adj = adjacency _matrix(tree) 
I = invariant-set^P) 
ol = tree_onebox(/, J 4c(7') 
achieved_isolation = (ol C S) 
S = oI 

tree = insert_onebox(tree,I) 

tree = insert_image(tree, f,oI) 

B = tree.boxnums() \ B 
until (achieved-isolation and B == 0) 
Output: S 



We can think of Algorithm 3 as a modification of grow isolating (Algorithm 1), 
where we make use of caching and lazy evaluation, meaning we only add boxes to 
the tree and compute box images when we absolutely have to. This way we can 
start with only our initial skeleton of points in the tree and grow both the tree 
and multivalued map just enough to accommodate the isolating neighborhood and 
verify its isolation. 

In fact we can see that we compute exactly the images we need, in a pre- 
cise sense. Consider the returned set S which, a-posteriori, must be equal to 
tree_onebox(J,v4(ij) for some invariant set /. Since we terminated, S must be the 
true one-box neighborhood of /; otherwise, we would have added new boxes in the 
insert_onebox call and failed to terminate. Similarly, the image of S — o(I) must 
already be in the tree, since otherwise the insert-image call would have resulted in 
new boxes. Thus, S is an isolating neighborhood of /, and since we grew S only by 
adding boxes to it, we have computed images only for boxes in S, which is precisely 
what we need to verify its isolation. 

Algorithm 3 is the cornerstone of our approach in this paper. Without it, no 
matter how clever our numerical approximations are, producing index pairs would 
essentially be just as expensive as computing the map on all boxes. This is es- 
pecially important when one is working with a low-dimensional invariant set in a 
high-dimensional embedded space, or when the hyperbolic, invariant sets whose ex- 
istence we wish to prove are tightly squeezed between non-hyperoblic sets or even 
singularities of the map. 

To see this more precisely, consider the memory required to store the box images 
using our bottom-up insertion approach as compared to the top-down approach of 
[DFT08] ; let Mj ns and Mdft be the memory in each case. As above let n be the 
dimension of the manifold, and K d be the number of boxes needed to cover the 
invariant set at depth d (which is independent of the method used). In settings 
where isolation is difficult, the method of [DFT08] would typically require -Mdft = 
0{2 dn ). Using Algorithms 2 and 3, however, we can achieve M ins = K d . Thus, 
in situations where K d -C 2 dn , we get a tremendous savings using our bottom-up 
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Figure 1. Plot of trajectories of the standard map for e — |. 

method, and if the depth required for isolation is high, this savings could be the 
difference between infeasible and feasible. We will see a concrete example of this 
in the next section when we study the standard map, with further discussion in 
section 4.4. 



4. Computations 

We apply our methods to the standard map 

/ e : (x, y) yx + y + — sin(27rx) mod 1, y + — sin(27r:r) mod 1^ , 

where e > is a perturbation parameter. This map is perhaps the best-known ex- 
act, area-preserving symplectic twist map. For e = every circle {y — constant} is 
invariant and the dynamics of the map consist rotations of this circle with frequency 

y- 

The map f e can be recovered by its generating function h : T 2 — > E in the sense 
that, if f e (x,y) — (X,Y), then YdX — ydx — dh(x 7 X). Moreover, we have its 
action A which, for a sequence of points {xjy, . . . , xm}, is 

Af-i 

(1) A{x N ,...,x M )= h(x k ,x k +i), 

k=N 

such that trajectories of f e "minimize" the action for fixed endpoints {x^, xm}, in 
an analogous way to classical Lagrangian mechanics (see [GolOl, §2.5]). 

It follows from a simple symmetry argument that for all e ^ 0, the stable and 
unstable manifolds of the hyperbolic fixed point intersect at x — \ (see Figure 2). 
Denote by (|,y e ) such a homoclinic point. It is rather easy to approximate this 
point numerically as one needs only to follow an approximation of either the stable 
or unstable manifold until it crosses x = \. Moreover, we can see from Figure 2 
that between the point (\,y e ) and its image f e (h,y e ) = (s + y s ,y e ) there is ho- 
moclinic point which we denote as (u,v). It is also relatively easy approximating 
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Figure 2 . Intersection of the stable and unstable manifolds at x = | . 



numerically this point based on having a good approximation of the stable and 
unstable manifolds and the point (^,y e ). Using these points, one can proceed in 
two different ways in order to get index pairs for and isolated invariant set which 
belong to the homoclinic tangle of the fixed point. 

The first way is to find the intersection of the stable and unstable manifolds at 
x = \ at high accuracy, and iterate this point enough times forward and backwards 
to get close enough to the fixed point. Considering all of these iterates and the fixed 
point as our invariant set, we may grow an isolating neighborhood and index pairs 
associated to them using Algorithm 3. This is perhaps the quickest way to get an 
index pair, as we already know where the isolated invariant set is. The downside 
is that in order for this approach to be successful, we need to compute the first 
homoclinic intersection with very high accuracy, as the homoclinic connection may 
be lengthy, and after enough iterations the error may become large enough to not 
give us a good approximation of the invariant set. 

The second approach is to compute the connections using shortest path algo- 
rithms as in Algorithn 2. This approach requires less precision in the computation 
of the homoclinic intersection, but is slower than the first approach. However, it 
is faster than a blind shortest-path search because it takes in consideration the 
dynamics of f e : recalling the action (1) of / E , we have A(xi,x 2 ) — h(xi,X2) and 
so we can compute the averaged action A : Q x Q — > R from box Bi to Bj as the 
average of A(x,y) with x £ \Bi\ and y £ \Bj\. Thus we can re- weigh the graph 
representing the map on boxes, from having every edge of weight 1, to having the 
edge going from Bi to Bj weight K + A(Bi,Bj), where K is any positive number 
satisfying 



K > K* = max \A(x\,x<2 
(xi,x 2 )eT 2 
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uniform for all Bi. Different choices of K do not necessarily give the same shortest 
(or cheapest) path: higher K gives more weight to the number of edges in the 
path (which one might use when working at a low depth), while lower K gives 
more weight to the action (which one might use at a high depth). Thus, searching 
for shortest paths in the graph with the new weights, we have a better chance to 
compute the right connection at the first try, as the connections which are computed 
contain, on average, the least action from the beginning box to the end box. 

This second method turns out to be a better fit in our setting, as numerically 
iterating the point (i,y s ) introduces error which is multiplied upon each iteration. 
Thus, the error grows exponentially fast until the iterates reach a small neighbor- 
hood of the fixed point, and this error is further exacerbated by the fact that we 
will mostly treat e as an interval. The second method still requires some precision, 
but it is modest enough that we can efficiently obtain it using the parameterization 
method [CFdlLOS]. Algorithm 4 summarizes our construction of index pairs which 
contain the homoclinic tangle of the hyperbolic fixed point. 

As e decreases, both the area of the Birkhoff zone of instability and the angle of 
intersection of stable and unstable manifolds of the hyperbolic fixed point (\,y e ) 
decrease exponentially fast with e [Gel9! I] . Thus since the intersection of the in- 
variant manifolds is barely transversal, the index pairs associated to the homoclinic 
orbits have to stretch out considerably across the unstable manifolds in order to 
achieve isolation (see Figure 4). This in turn implies that Algorithm 3 takes more 
iterations to cover the invariant set. Moreover, the size of the boxes is necessarily 
exponentially small with e, and so the number of boxes needed to create the index 
pairs increases rapidly. The bottom line is that the complete automation of the 
procedure prevents us of having to create the index pairs by hand, which for low e 
must be an extremely difficult task. 

KAM theory asserts that for |e| small enough (roughly |e| < .971 [Jun91]), there 
is a positive measure set of homotopically non-trivial invariant circles on which 
the dynamics of f e is conjugate to irrational rotations. In this case the invariant 
circles folliate the cylinder and serve as obstructions to orbits from wandering all 
over the cylinder, i.e., each orbit is confined to an area bounded by KAM circles. 
In this case, the topological entropy of f £ is concentrated in the Birkhoff zone of 
instability associated to the homoclinic tangle of the hyperbolic fixed point. Once 
e > e* « .971 there remain no homotopically non-trivial KAM circles to bound the 
y coordinate of the orbit of a point and one has then hope to find connecting orbits 
between different hyperbolic periodic orbits. 

We apply our methods to obtain three types of results: 



• For all f E with e G [-7,2], we give a positive lower bound for its topolog- 
ical entropy. This is done by treating e as an interval. An advantage of 
having the procedure automated is that one can easily study a parameter- 
depending system at different values of the parameter. Treating the pa- 
rameter as an interval allows us to detect behavior which is common to all 
values of the parameter in the interval. This is done in section 4.1. 

• Not treating e as an interval allows our method to go further and obtains 
positive bounds for the much lower value of e = \. This is done in sec- 
tion 4.2. 
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• Our examples in section 4.3 combine the new methods of this paper with 
the spirit of [DFT08] of connecting periodic orbits to find better entropy 
bounds, which will illustrate for the case e = 2. 

We remark that in [MJ09, §4.1] an alternate approach for the creation of index 
pairs for the standard map is given. It is done through set-oriented methods which 
are based on following the discretized dynamics along the discretized stable and un- 
stable manifolds. We do not know how this approach would perform when treating 
e as an interval, although we suspect it would perform equally well. Our bottom-up 
approach for constructing index pairs requires the computation of fewer box im- 
ages, but in general we may make use of slightly more a-priori information such 
as the knowledge of where hyperbolic invariant sets are. The algorithms in [MJ09] 
(and indeed those of [DFT08, §3.1]) require less a-priori knowledge of hyperbolic, 
invariant sets, but require a greater number of box-image computations (see [MJ09] 
for more details). 

We should point out that the bounds we provide in the following sections are 
close to some of the non-rigorous bounds given in [YT09]. To our knowledge there 
are no other bounds in the literature for the topological entropy of the standard 
map for small values of e. It is expected that the entropy is exponentially small 
as e — > [Gel99] , while it is known that the entropy grows at least logarithmically 
in £ as e — > oo [Ke ]. But for small values of e, we have not been able to find 
computational bounds besides the ones already cited. 

4.1. Parameter exploration. As remarked earlier, when e < e* there exist in- 
variant KAM circles which prevent the connections between many periodic orbits. 
This forces us to concentrate on the homoclinic tangle of the hyperbolic fixed point. 
The results from this section are obtained using Algorithm 4 and summarized in 
Theorem 4.1. 



Algorithm 4 Creating index pairs for f e using the homoclinic orbits of the fixed 
point 

Input: e = [e~,£ + ], fg 

Let yg = and compute (5,2/e) using [CFdlL05]. 

H = {(\,Ve)\ (f , 1 " Ve) , (u, »), (1 - u, 1 - v), (0, 0)} 
X = [j p ^ qeH {p,q) 

Find connections using Algorithm 2 and the weighted graph using the averaged 
action A. 

Grow the index pairs using Algorithm 3. 



By performing the computations using e as an interval, we are proving behavior 
which is common for all f e within such interval. In such case then our guess for 
the homoclinic intersection (\,y e ) is done only for one point in the interval (the 
midpoint). In general it is easier to isolate an invariant set for smaller parameter 
intervals, since the wider the interval, the more general the isolation must be. In 
our setting, it is much easier to isolate homoclinic connections for e > e* than it 
is for s < e*. To reflect this, using a crude approximation of e* ~ 1.0 for ease of 
bookkeeping, we use intervals of size 0.005 for e > 1.0, but we shrink our interval 
size to 0.001 for e < 1.0. 
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As e decreases, the size of the boxes we use decreases, and the length of the 
homoclinic excursion increases, leading to another increase in the number of boxes 
needed and a longer running time of Algorithm 3. At some point, it becomes 
computationaly unrealistic to continue; we stopped somewhat before this point, 
when the interval computations took roughly 40 hours for e = [0.700,0.701]. See 
section 4.4 for further discussion of the implementation and efficiency. 

Theorem 4.1. The topological entropy of the standard map f e for e € [-7,2] is 
bounded from below by the step function given in Figure 3. In particular, we have 
h(f £ ) > 0.2 for all e £ [.7, 2]. The precise individual values for each subinterval are 
given in Appendix A. 
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Figure 3. Lower bounds of h(f e ) as a function of e in the interval 
[.7,2]. Note that all bounds in this interval exceed 0.2. 

Proof. For each of the e-intervals e = [e~ , e + ] on the table found in Appendix A, we 
get an index pair for an isolated invariant set of fg using Algorithm 4. Using then Al- 
gorithm 5, 6, 7, and 8 and Theorem 3.6 from [DFT08], which is essentially amounts 
to a finite number of checks using Corollary 2.13, we prove a semi-conjugacy to a 
subshift of finite type, from which we get a bound on the entropy by bounding the 
spectral radius of the associated matrix. □ 

There are three aparent scales on which our lower bounds for h(f £ ) change with 
respect to e: global, semi-local, and local. Clearly there is an evident global increase 
of the entropy bounds as e increases, as is to expected. On a semi-local level, there 
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are a few intervals (roughly [1.6, 2], [1.2, 1.6], [.92,1.2], [.8, .92], [.75, .8]) on which 
the bounds for h(f e ) seem to hover around a fixed value per interval. This is due 
to using the same depth on such intervals. As e decreases, we need to increase the 
depth. Locally, the apparent irregularity of the function of lower bounds is due to 
the nature of our automated approach: the accuracy of the guess ( \ , yg) varies per 
interval, as does the computation of the averaged action, et cetera. 

4.2. Positive bound of h(f £ ) for lowest e. In this section we show an example 
of an index pair for the standard map for e = A . This value was picked because it is 
small enough that we can illustrate the strengths of our algorithms in tight places 
while keeping the computation times reasonable. Using Algorithm 4 with e = | we 
get an index pair shown in Figure 4 (although it is barely visible). 

Theorem 4.2. The topological entropy for the standard map f e when e = \ is 
bounded below by 0.1732515918346. 

The proof is the same as in Theorem 4.1. The tree from which the index pair 
obtained for Theorem 4.2 was obtained contains 568,754 boxes. Among those, 
281,530 are in the index pair. Roughly a quarter of the boxes in the index pair 
form the exit set (64,518). This index pair (Pi,P ) gives us an induced map which 
acts on Hi (Pi, P ; Z) = Z 1801 but is reduced to a SFT in 73 symbols. 

Figure 4 shows, on top, the plot of some trajectories for the standard map at 
e = |, of the stable and unstable manifolds, and an index pair for the homoclinic 
orbits. On the bottom is a close-up to a component of the index pair squeezed by 
KAM tori and a barely transversal intersection of the stable and unstable manifolds. 
The boxes making up this index pair are of sides of size 2 -15 . The strength of our 
"growing-out" approach is that the creation of such index pairs in tight places is 
achievable and can be automated. 

4.3. Higher periods. When e > e* all homotopically non-trivial KAM circles are 
vanished, thus it is possible to connect different periodic orbits. The appendix 
of [Gre79] contains an algorithm for finding periodic orbits for the standard map. 
We implement this method to find the periodic orbits which we use to grow index 
pairs. 



Algorithm 5 Creating index pairs using the homoclinic orbits of the fixed point 
and periodic orbits 

Input: e= [£",£+], f e , PeN 

Let y§ = y s - +s + and compute (\,yg) using [CFdlL05]. 
Hi = {(0, 0), 2 (i, fe) , (|, 1 - ife) , («, v), (1 - u, 1 - «)} 

H2 = hyperbolic periodic orbits of fg up to period P (computed using the ap- 
pendix in [Grc79]) 

Find connections using Algorithm 2 and the weighted graph using averaged action 
A. 

Grow the index pairs using Algorithm 3. 



Algorithm 5 is essentially the main strategy employed in [DFT08]. In that paper, 
good index pairs were found by finding pairwise connections between periodic orbits. 
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Figure 4. On top, a plot of the trajectories of the standard map for e — ^ along 
with the stable and unstable manifolds of the hyperbolic fixed point which seem to overlap. 
On the bottom, a close-up of a component of the index pair yielding Theorem 4.2 which 
contains the homoclinic point at x — i along with trajectories, some of which are part 
of KAM circles squeezing the component. It overlays the stable and unstable manifolds, 
whose angle of intersection is very small, causing the index pair to be very sheared. 

Such an approach can be slightly generalized by looking for pairwise connections 
between hyperbolic, invariant sets, which is what Algorithm 5 does. 

We apply the algorithm to e = 2.0, with maximum period P — 2. Figure 5 
shows the index pair for this computation. We remark that besides finding pair- 
wise connections between hyperbolic periodic orbits, we find connections between 
periodic orbits and the homoclinic orbit (|, y e ) and (w, v) mentioned in section 4.1. 
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This allows us to find richer dynamics and to achieve higher entropy bounds. The 
result from this index pair is summarized in the following theorem. 

Theorem 4.3. The topological entropy for the standard map f e fore = 2 is bounded 
below by 0.44722970117798. 

The proof is again similar to Theorem 4.1. The index pair (Pi,Pq) has a total 
of 8600 boxes, and gives us an induced map which acts on H 1 (Pi, P ;Z) = Z 105 
which is reduced to a SFT in 59 symbols. 



I* 



w 



WW 



09 1+ I 4 

0.8 - Y*A^ "^^^ 

07- \ 



.i 



0.1 -V ^ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.S 0.9 1 

Figure 5. Index pair obtained through Algorithm 5 for e = 2 and 
P = 2 at depth 9. 

This result shows that our bounds in Theorem 4.1 are probably suboptimal, 
since our interval bound for [1.995, 2.000] is lower than in Theorem 4.3. It is not 
surprising that by adding connections to another hyperbolic periodic orbit we may 
find a higher entropy bound. 

Algorithm 5 can be quite effective when P is small, but as P grows, there is large 
increase (quadratic at the very least, often exponential) in the number of connec- 
tions to compute, especially since the number of period- P orbits grows rapidly with 
P for chaotic maps. Moreover, the number of long connections increases, which as 
discussed in section 3.2 makes Algorithm 2 work even harder. As suggested in that 
section, we instead turn to computing periodic orbits of higher period rather than 
computing connections explicitly. 

For our last example, we apply this approach on top of the previous example 
from Theorem 4.3 to produce a very strong index pair. That is, we run Algorithm 5 
with P — 2 and then simply add all hyperbolic orbits of period 3 and 4 (there are 
four of each) without adding any further connections. The resulting index pair, 
shown in Figure 6, clearly benefits from these "natural" connections. 

Theorem 4.4. The topological entropy for the standard map f £ fore — 2 is bounded 
below by 0.54518888942276. 
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Figure 6. Index pair obtained by adding periodic orbits of period 
3 and 4 to the example in Figure 5 

The index pair (Pi,Po) for Theorem 4.4 has a total of 64,185 boxes, with 5,839 
in the exit set. The induced map acts on H\{P\,Pq\'L) = Z 138 and reduces to a 
SFT on only 41 symbols. 

4.4. Notes on efficiency and implementation. The computations in this paper 
were performed in MATLAB on machines with between 1 and 2 gigabytes (GB) 
of memory and with clock speeds between 1.5 and 2.5 gigahertz. Runtimes ranged 
from 3 or 4 minutes for e = [1.995,2] and e — 2.0, to almost 2 days for e = 
[0.700,0.701] and roughly 5 days for e — 0.5. As discussed in section 4.1 and the 
beginning of section 4, there was a roughly inverse exponential relationship between 
e and the runtime for the intermediate values. 

The two most time-consuming subroutines for our computations were Algo- 
rithms 2 and 3, which of course is one reason for our focus on them in this work. 
As e decreased, however, Algorithm 3 dominated the runtime, particularly in the 
bookkeeping step (maintaining the correct box numbers, as discussed in section 3.1) 
and the insertion step, when new boxes are inserted into the tree. While the inser- 
tions cannot be avoided, this does suggest that a better tree implementation could 
enhance performance greatly. 

To conclude, we would like to reiterate how difficult it would be to reproduce our 
results for low e, namely Theorem 4.2 and the lower intervals of Theorem 4.1, using 
index pair algorithms from [DFT08]. As mentioned in section 3.3, our algorithms 
are considerably more memory efficient in certain situations, which we now explore 
concretely. Consider the index pair we obtained in Theorem 4.2 when e = 0.5; 
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there were 568754 boxes in the index pair, at depth 15, and the adjacency and 
transition matrices took up about 0.1GB of memory using our approach. Using 
the [DFT08] approach, one would need to compute the map on all boxes at depth 
15. Using a conservative estimate of 190 bytes per box to store the image and 
adjacency information (the average for our index pair was 194.38 bytes), this would 
require 204GB, which is beyond reasonable at the time of this writing. More to 
the point, our bottom- up approach is clearly orders of magnitude more efficient in 
terms of memory. When one considers the graph computations that would need to 
be carried out on the resulting 2 30 -node graph (the transition matrix), it becomes 
even clearer that our computations would have been impractical using the top-down 
approach of [DFT08]. 



Appendix A. Precise bounds for h(f e ) 



Below we list the actual values for the lower bounds in Theorem 4.1 the last 
column indicates the number of symbols for each SFT whose entropy bounds that 
of f E for e in each interval. 
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